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The promise of single Cooper pair quantum circuits based on tunnel junctions 
for metrology and quantum information applications is severely limited by the 
influence of "offset" charges - random, slowly drifting microscopic charges 
inherent to many solid-state systems. By shunting a small junction with the 
Josephson kinetic inductance of a series array of large capacitance tunnel 
junctions, thereby ensuring that all superconducting islands are connected to 
the circuit by at least one large junction, we have realized a new superconduct- 
ing artificial atom which is totally insensitive to offset charges. Yet, its energy 
levels manifest the anharmonic structure associated with single Cooper pair 
effects, a useful component for solid state quantum computation. 

Electric charge can be manipulated at the level of a single charge quantum (1) in two types 
of superconducting circuits with different topologies. The minimal example of the first type of 
circuit is the Cooper pair box, which consists of an isolated superconducting electrode ("island") 
connected to a superconducting reservoir on one side by a small tunnel junction, and on the 
other side by a gate capacitance in series with a voltage source. The dynamics of the island is 
described by two variables: the integer number of Cooper pairs occupying the island and its 
conjugate, the 27r-cyclic superconducting phase difference between the island and the reservoir. 
The junction area must be chosen sufficiently small such that the electrostatic energy of the 
island due to an extra Cooper pair is larger than the Josephson energy of its coupling to the 
reservoir, thus confining fluctuations of the number of Cooper pairs below unity. Stated in 
electrical engineering language, one needs Zj > Rq, where the junction reactive impedance 
Zj = {Lj/CjY^"^ is defined by the Josephson characteristic inductance Lj and capacitance 
Cj (2), and where the superconducting impedance quantum is given by Rq = h/ (2e)^ ~ 1 kfi, 
denoting Planck's constant h and the charge quantum e. The second type of circuit is based on 
a superconducting loop connecting the two electrodes of a small junction with an inductance 
which exceeds Lj. The circuit conjugate variables are now the magnetic flux generated by the 
persistent current in the loop and the displacement charge on the plates of the small junction 
capacitance. When Zj > Rq, the large loop inductance is submitted to quantum fluctuations 
of flux larger than the flux quantum $o = 2nh/2e, and therefore according to Heisenberg 
principle, the junction charge fluctuations are reduced below the value 2e. 

In practice, the realization of both circuit types faces fundamental difficulties. Islands are 
exposed to random electric fields due to fluctuating charged impurities which are ubiquitous 
in most solid-state environments and whose compounded effect is described by a noisy offset 
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charge. Although the fully developed charging effects were demonstrated for the Cooper pair 
box (3, 4), it soon became clear that the low-frequency offset charge noise was a major source of 
decoherence for charge qubits derived from this device (4-7). This state of affairs has prompted 
the development of alternative superconducting qubits based on large junctions with Zj <^ Rq, 
avoiding the single Cooper pair regime and the related charge offset problem (8-10). On the 
other hand, implementing the island-free circuit, which is immune to charge offset noise, is an- 
other hard problem. This is because any finite-length wire with inductance L always comes with 
self -capacitance C which reduces the total charging energy of the circuit and therefore steers 
it away from the charging regime, unless (L/C)^/^ ^ Rq. In fact, a purely electromagnetic 
inductance is incompatible with the single Cooper pair effects since (L/C)^/^ is then bounded 
by the vacuum impedance (/xo/^^o)^^^ ~ 377 H < Rq, hq and Eq being vacuum permeability 
and permittivity (11, 12). 

In this paper, we present experimental results on a novel single Cooper pair circuit based on 
a superconducting loop, which solves both the inductance and the offset charge noise problems. 

The small junction of our circuit is shunted by a series array of carefully chosen larger 
area tunnel junctions (Fig lA-C). Here, all islands are connected to the rest of the circuit by 
at least one large junction so that quasistatic offset charges on all islands are screened. The 
large capacitances of the array junctions prevent phase slips within the array, and for excitations 
whose frequencies are below the junction plasma frequency, the array effectively behaves as 
an inductive wire. By choosing a sufficiently large number of array junctions it is possible to 
create an inductance exceeding that of the small junction. At low energies, the loop is effectively 
described by the loop flux $ and the small junction charge Q, satisfying [!>, Q] — ih. 

To form a charge offset-free inductively shunted junction, four conditions involving the ef- 
fective inductance Lja and capacitance Cja of the Af array junctions are required: (i) MLja ^ 
Lj, (ii) e-^^'Ql^-'^ < e < 1, (iii) Ue-^^'-^'^-'^ < e-^^-Q'^J, and (iv) M < {CjAlCgfl^. In 
the first relation (i), we simply estimate the total array inductance to be MLja and require that 
it exceeds the small junction inductance, allowing it to support the large flux fluctuations of 
the loop. The second relation (ii), where Zja — (Lja/Cja)^/'^ is the array junction reactive 
impedance, dictates the minimum size of the array junctions necessary to reduce (13) the un- 
controlled offset charge on the islands of the circuit below the desired value of the order of 
2e X e. The third relation (iii) ensures that the inductive role of the array is not jeopardized by 
quantum phase slips (14). Specifically, the probability amplitude of a phase slip event within 
the array (l.h.s.) must be negligible compared to that in the small junction (r.h.s.). According to 
relation (iii) a fluxon tunnels in and out of the loop predominantly via the small junction, thus 
effectively erasing the discrete character of the array. Lastly, relation (iv) states that the induc- 
tance of the array is not shunted by the parasitic capacitances Cg of array islands to ground. It 
is obtained by estimating the array parasitic resonance frequency to be {LjaM x CgAf)~^^'^, 
and requiring that it is larger than the junction plasma frequency (LjaCja)^^^"^- Remarkably, it 
is the relation (iv) which, with present junction technology, most severely limits the maximum 
number of junction in the array and, thus, its maximum inductance. 

We have implemented the above array proposal and constructed a new superconducting arti- 
ficial atom which we have nicknamed "fluxonium". It contains M — A3 Al-AlOx-Al Josephson 
junctions (15) such that Zja 0.5Rq and a small junction with Zj ~ 1.5Rq (16). The above 
four conditions being realized, the fluxonium can be modelled (Fig. ID) as a small junction 
shunted by an inductance La (17). The three characteristic energies of this model, namely 
El = (^q/2tiY /La, Ej = ($o/27r)^/Lj and Eq = e^/ (2Cj) have values corresponding to 
0.52 GHz, 9.0 GHz and 2.5 GHz, respectively. The additional L^Cr resonator, capacitively 



2 




Figure 1: (A) Sketch of a small Josephson junction shunted by an array of larger area junctions. 
The two superconducting leads of the small junction are coupled capacitively to a quarter-wave 
microwave resonator, a "parallel wire" transmission line shorted on the opposite end. The 
resonator itself is probed capacitively and symmetrically via two 50 microwave ports, result- 
ing in a loaded quality factor of 400. The whole device is made with a single step standard 
Al/AlOx/Al double angle evaporation through an e-beam lithography mask on a high resistivity 
Si substrate. (B) Close-up view of the small junction region, showing top and bottom junction 
electrodes (grey) and their thin oxide layer (green). Array junctions are about one order of mag- 
nitude larger in area and spaced as tight as e-beam lithography resolution allows, minimizing 
microwave parasitics. (C) Electrical circuit representation of the loop formed by the small junc- 
tion (black), with Josephson inductance L j and capacitance Cj, shunted by the array of larger 
junctions (purple), with the corresponding inductance Lja and capacitance Cja- Islands formed 
between the array junctions have small capacitance to ground Cg. (D) Simplified circuit model 
of the fluxonium, consisting of three sections: i) circuit equivalent to a Cooper pair box, where 
the small junction with capacitance Cj and non-linear Josephson inductance Lj is capacitively, 
with capacitance Cc, coupled to the probe (solid black), such that Lj/CjY^"^ > h/{2eY ; ii) 
giant inductance La ^ Lj provided by the junction array (purple); iii) a parallel combination 
of Cn and such that {L^/CjiY^'^ ~ 50i7 <^ h/(2eY which is the circuit model for dis- 
tributed transmission line resonator (grey). (E) Potential seen by the reduced flux ip and energy 
spectrum of the circuit (D) for two values of external flux $ext- At $ext = 0, energy levels pos- 
sess well-defined parity as indicated with and '-' signs next to level numbers. Note that, in 
contrast with the RF-SQUID or flux qubit, there is on average only one level per local minima. 
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Figure 2: Modulation of the reflected 8.18 GHz microwave signal with externally applied flux 
$cxt- The signal is clearly flux-periodic indicating that the junction ring is closed and super- 
conducting. The values of $cxt at which the signal undergoes full swings correspond to the 
anticrossings of the — 1 transition frequency of the device with the resonator bare frequency, 
later inferred to be 8.1755 GHz. The measurement tone populates the resonator with less than 
0.01 photon on average. 

connected to the small junction (Fig. ID), reads out the "atom" in a manner analogous to the dis- 
persive measurement of cQED qubits (18). It is implemented by a quarter-wave superconduct- 
ing coupled microstrip resonator (Fig. lA) with quality factor 400 due to capacitive coupling to 
the two 50 f2 measurement ports. The resonator frequency cor = {LjiCji)^^^'^ ~ 27r x 8.17 GHz 
is pulled by the reactance of the fluxonium circuit and is monitored by standard ultra low noise 
microwave reflection technique. The fluxonium reactance depends on its quantum state, an 
effect leading to a purely dispersive state measurement (15). An externally imposed, static 
magnetic flux $cxt threading the loop •I'o-periodically modulates the spacings of energy levels 
of our artificial atom. 

Introducing the operators N = Q/2e and (p = 2e^/h, describing the reduced charge on 
the junction capacitance and its conjugate reduced flux operator (i), the Hamiltonian of the 
fluxonium coupled to its readout resonator can be written as 

H = AEcN^ + ^^^(^2 _ _ 27r<l>e,t/$o) +gN{a + a^) + hwRa^a (1) 

Here a is the photon annihilation operator for the resonator, g is the atom-resonator coupling 
constant. The second term and the range of definition of and N, whose eigenvalues are here 
both on the entire real axis, distinguishes the form of Hamiltonian (1) from that of the Cooper 
pair box in cQED experiments {18). There are three important points to note concerning this 
Hamiltonian (2): i) it is invariant under the transformation N ^ N + A^oflfset (A^offset stands 
for offset charge value) hence the "charge-free" character of our device; ii) it differs from that 
of the transmon {13) since offset charge influence is screened for all states, not just for the 
low-lying states; iii) its second term, despite the fact that is the smallest of the fluxonium 
energies, has a non-perturbative influence on the full energy spectrum of this artificial atom, 
which presents strongly anharmonic transitions {21) (Fig. IE). Our experiment probes these 
transitions by microwave spectroscopy, from which we infer the size of charge fluctuations. 

To characterize the fluxonium, we first measure the ground state resonator pull as a function 
of <l>ext- The results (Fig. 2) show the expected $o~ periodicity as well as the avoided crossings 
of the resonator frequency and the ground to excited state transitions. This confirms that the 
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Figure 3: Phase (colored circles) of reflected readout tone as a function of spectroscopy tone 
frequency taken at $cxt = 0.05$o- Data for the first three resonances (further identified as 
transitions from the ground state to states 1, 2 and 3) is shown from left to right in red, orange 
and green respectively. Resonances are well-fitted by Lorentzians (solid black lines) for a broad 
range of spectroscopy powers. Insets on the two sides of each resonance show the dependence 
of the resonant peak height (left) and width squared (right) on the spectroscopy tone power. Data 
in all insets follow the predictions (solid black lines) of Bloch equations describing relaxation 
dynamics for a spin 1/2 and indicate that all transitions involve one photon. 

entire 44 junction loop is superconducting and that the resonator-atom system is in the strong 
coupling regime of cavity QED (22). 

Next, we perform a two-tone spectroscopy measurement (2i) at a fixed flux $cxt = 0.05$0' 
during which, in addition to the fixed frequency readout tone, we probe the transition frequen- 
cies of the atom through a second, variable frequency spectroscopy tone. The resulting peaks 
(Fig. 3), correspond to the later-determined — 1,0 — 2, and — 3 transitions from the atom 
ground state. The peaks are well-fitted by Lorentzians and their power-dependent widths and 
heights are well-explained by the Bloch equations of precessing spin 1/2 (24) as shown in the 
insets of (Fig. 3). Extrapolating fitted linewidths to zero spectroscopy power, we obtain lower 
bound estimates of their decoherence time at 350, 250 and 80 ns respectively. 

Our main result is the spectroscopic data collected as a function of both spectroscopy fre- 
quency and flux (Fig. 4A). Note that $cxt variations span 20% of $o around $cxt = instead 
of the usual 1% or less around $o/2 in flux qubit experiments (9). In (Fig. 4B) we compare the 
measured peak center frequencies with the prediction for the — 1,0 — 2, — 3 and the two- 
photon — 4 transitions obtained from numerical diagonalization of Hamiltonian (Eq. 2). Note 
that we are in effect fitting more than three flux dependent functions, i.e. the flux dependent 
transition frequencies, with only three a priori unknown energies Eq, El and Ej so the problem 
is severely overconstrained. The fit of the line (Fig. 4B) labeled SR (for array self -resonance) 
requires a minor extension of the model taking into account parasitic capacitances across the 
array (15). Apart from introducing another resonator mode coupled to the atom, this extension 
by no means invalidates the inductive character of the array, at least as far as the — 1 and — 2 
transition of the fluxonium are concerned. Even the perturbation of the — 3 and — 4 transition 
frequencies by this extra mode is less than 2%. 

Based on the excellent agreement between theory and experiment, we infer the wavefunc- 
tions of the first three energy levels, and plot their amplitudes both in charge (Fig. 4C) and flux 
(Fig. 4D) representations for $ext = 0. In the ground state, we find that the ratio of charge 
to flux fluctuations is AN/ Acp = 0.56, about 5 times smaller than the fine structure constant 
allows for a conventional resonator. This confirms that the charge in our circuit is indeed lo- 
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Figure 4: (A) Phase of reflected readout tone as a function of the spectroscopy tone frequency 
and external flux. The color scale encodes the value of phase, zero corresponding to the mauve 
background, blue to positive values (peaks) and red to negative values (dips). The grey re- 
gion shows the reflected phase of a single tone, swept close to the resonator bare frequency 
exhibiting a 50 MHz vacuum Rabi splitting of the resonator with the fluxonium transition — 1. 
The inset in (A) zooms in the central region of the — 2 transition and confirms that it is in- 
deed symmetry-forbidden at $cxt = 0. (B) Measured peak frequencies (blue circles) fitted 
by numerically computed spectrum of Hamiltonian (2) (solid red lines) and its modification 
(see supplementary online text) to explain the additional transition labeled "SR" (dashed black 
lines). (C) Amplitude of fluxonium wave functions for levels (black), 1 (red) and 2 (orange) 
computed in charge representation at zero flux bias, using circuit parameters extracted from 
the fits. (D) Same as in (C) but in flux representation. The flux representation wave func- 
tions demonstrate that the reduced flux is delocalized compared to the size of the Josephson 
well while charge wave functions confirm that the localization of charge on the junction is less 
than a single Cooper pair charge. Note that in this circuit, the junction charge is a continuous 
variable, in contrast to the Cooper pair box, and flux swings of more than 27r are allowed. 
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calized at the single Cooper pair level (AN — 0.53, Aip = 0.95). The wavefunctions in flux 
representation (Fig. 4D) can be interpreted as simple superpositions of states in which the re- 
duced flux if is localized in the wells of the Josephson cosine potential (fluxon states, hence the 
name fluxonium). The parity of fluxonium states, which forbids the — 2 transition at zero ex- 
ternal flux, manifests itself explicitly by a remarkable "hole" in the corresponding spectroscopic 
line (Fig. 4A, inset). The allowed transition between the second and third level (data not shown) 
is particularly spectacular since it corresponds to motion of the total flux in the fluxonium loop 
by two whole flux quanta. This is to be contrasted with the 10% of flux quantum or less flux 
motion involved in transitions of the flux and phase qubits (8, 9). Nevertheless, despite the large 
flux fluctuations of the system and the corresponding charge pinning, the circuit has complete 
immunity to offset charge variations: the data of (Fig. 4A) has been taken piecemeal in 72 hours 
and no jumps or drifts have been observed during this period. 

We have thus demonstrated that an array of Josephson junctions with appropriately chosen 
parameters can perform two functions simultaneously: short-circuit the offset charge variations 
of a small junction and protect the strong non-linearity of its Josephson inductance from quan- 
tum fluctuations. The data shows that the array possesses a microwave inductance 10^ times 
larger than the geometric inductance of a wire of the same 20 //m length. The reactance of such 
inductor is about 3Rq ~ 20 kQ at 10 GHz while its resistance is less than 1 il. The spectrum 
of the fluxonium qubit suggests it is as anharmonic as the flux qubit but as insensitive to flux 
variations as the transmon qubit. Possible applications of this single Cooper pair charging effect 
immune to charge noise include the observation of fully developed macroscopic quantum coher- 
ent oscillations between fluxon states (25), the search for a "A" or "V" transition configurations 
for the shelving of quantum information (26) in superconducting artificial atoms, topological 
protection of superconducting qubits (27), and finally the long-sought quantum metrology of 
electrical current via Bloch oscillations (28, 29). 
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1 Materials and Methods 

Sample fabrication. The device is made on a high-resistivity Si substrate, 300 jim. thick. Both 
Josephson junctions and the readout resonator are fabricated in a single step using e-beam 
lithography, double angle Al e-beam evaporation and lift-off techniques. The Al evaporation 
and oxidation is conducted in an e-gun evaporator at pressures less than 10^^ Pa, AlOx grown 
in the environment of 680 Pa of 15% oxigen-in-argon mixture during 10 minutes. The areas 
of the small junction and array junctions are designed to be 0.2 x 0.3 //.m^ and 0.25 x 2 /xm^, 
respectively. All 43 array junctions are equally spaced at less than 200 nm so that total length 
of the array is only 20 /xm. The loop area of the array- small junction ring is 3 x 20 jiw?. 

Sample mount. The Si chip is glued using GE varnish to the copper base of a fully enclosing, 
custom-made microwave sample holder, shielding the sample from both residual RF, infrared 
and optical photons. The holder provides two well-matched transitions from the Anritsu K- 
connectors on the outside of the holder to the two microstrip lines made on a PCB inside the 
holder. The resonator's on-chip launching pads, schematically indicated as two sections of a 
coaxial cable and marked "50 Q" in (Fig. lA), are then wirebonded to the ends of the two 
microstip lines. 

Cryogenic setup. The experiment is performed in a dilution refrigerator with base temper- 
ature 10 — 20 mK. Both resonator and the qubit are differentially excited via the A-port of a 
180° hybrid (Krytar, 2 — 18 GHz), whose two outputs are connected to the two ports of the 
sample holder. Incoming and outgoing signals are separated with a directional coupler (Krytar, 

2 — 20 GHz). The incoming signal line is attenuated using 10 and 20 dB microwave attenuators 
(XMA) at all temperature stages of the refrigerator, to remove non-equilibrium noise. The out- 
put line is amplified at the 4 K stage with a low-noise HEMT amplifier (Caltech, 1 — 12 GHz, 
30 dB gain). Two cryogenic isolators (Pamtech, 4 — 12 GHz, 15 dB) are placed between the 
amplifier and the sample, at the 800 mK stage and at the base stage, again to remove non- 
equilibrium noise, especially that coming from the amplifier. Stainless steel SMA cables are 
used to connect between the different temperature stages. All components are thermally an- 
chored to the proper refrigerator stages. A ~ 1 cm diameter custom made superconducting coil 
is glued to the sample holder, a few mm away from the chip, to provide perpendicular magnetic 
flux bias. The sample holder together with the coil is placed into a Cryoperm cylinder to shield 
it from stray quasistatic magnetic fields. 

Room temperature measurement setup. The readout resonator is excited with Agilent E8257D 
signal generator, the spectroscopy signal is generated using Agilent E8267D vector signal gen- 
erator and Tektronix 520 AWG. Both signals are combined at room temperature and sent into 
the input line of the refrigerator. The reflected ~ 8 GHz readout signal from the refrigerator out- 
put line is amplified at room temperature with two Miteq amplifiers (1 — 12 GHz, 30 dB gain), 
mixed down with a local oscillator (a third Agilent E8257D) to an IF signal of — 50 MHz, 
filtered and amplified with the IF amplifier (SRS SR445A), and finally digitized using 1 GS/s 
Agilent Acqiris digitizer. A software procedure then extracts the phase and the amplitude of the 
digitized wave. The experiment is typically repeated 10^ times to average the Gaussian noise to 
an acceptable level. Because the duration of each experiment is about 10 microseconds, every 
averaged data point is taken in a fraction of a second. All microwave test equipment is phase 
locked using a Rb precision 10 MHz reference (SRS FS725). The magnetic coil is biased in 
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Figure 5: Models for the fluxionium device. (A) Inductively shunted junction model. (B) 
Extended fluxionium model, including the capacitive coupling to the mode of a transmission- 
line resonator, and parasitic capacitances across the array. Numbers in cyan enumerate the 
nodes of the circuits. 

series with a resistor with Yokogawa 7751 voltage source. 

Comments on the data. The data in (Fig. 2) shows the digitized homodyne (zero IF) signal 
as a function of magnetic field, with the spectroscopy generator turned off. The data in (Fig. 
3) shows the phase of the digitized heterodyne (50 MHz IF) signal, as a function of frequency 
of the spectroscopy generator. The data in (Fig. 4A) is taken in the pulsed regime, when the 
spectroscopy generator outputs a 6 /Us saturating pulse followed immediately by the 2 fis readout 
pulse. This way we ensure that the sample is exposed to only one tone at a time, avoiding various 
spurious effects. The image presented in (Fig. 4A) contains 367 x 4597 data points. 



2 Supplementary Text 

In our analysis of the fluxonium device, we use two simple models whose corresponding cir- 
cuits are depicted in Figure [51 (A) the inductively shunted junction model, (B) the extended 
fluxonium model, describing the fluxonium coupled to a transmission-line resonator. 

Inductively shunted junction model. The simplest model of the fluxonium device is the 
inductively shunted junction model, see Fig.[5tA). It neglects parasitic capacitances across the 
fluxonium's Josephson junction array, and assumes that all internal degrees of freedom of the 
array are frozen out. In this limit, the fluxonium consists of a single Josephson junction with 
capacitance C j and Josephson energy E j, shunted by a large inductance L. Quantization of this 
circuit (1) is straightforward and leads to the Hamiltonian 

^0 = ^EcN^ + ]^ELf - Ej cos - 2vr^^ , (2) 

where the charge on the junction capacitance N (in units of 2e) and the reduced flux are 
canonically conjugate variables, [(^, N] = i. Structurally, this Hamiltonian is identical to the 
Hamiltonian describing one-junction flux qubits, and flux-biased phase qubits. However, the 
regime of large inductances relevant for the fluxionium differs from typical parameters in flux 
and phase qubits, and has been discussed in Ref. (2). 

Extended fluxonium model. For a more complete modelling of the spectra obtained in the 
experiment, we take into account the coupling of the fluxonium device to a transmission-line 
resonator. In addition to this resonant mode, the experimental data shows another resonance 
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coupling to the fluxonium device. Such additional resonances are expected when accounting 
for the parasitic capacitances of the Josephson junction array. The simplest effective model 
accurately describing the experimental data includes very few of these capacitances, and ap- 
proximates the array by a combination of inductances and capacitances as shown in Fig.[2B). 
The Lagrangian describing this circuit can be written in the form 

^ '^a - 2l(03 - $ext)^ + Ej COS + - i^ei + - sfcv'l (3) 

+ f 0304 + 

where C'j = CJ2 + Cj + 2CiXl + C^Xl Xi = U/L, C = 9(Ci + 2C2)XlXl, L' = 1/(9X1X2), 
C'ji = Cr + Cc/2, C = 6AiA2(CiAi — C2X2), and we have disposed of another resonant mode 
which does not couple to the fluxonium device. In terms of the original generalized flux (pi 
at each node i, the relevant variables are ^3 (associated with the fluxonium subsystem), the 
resonator mode (/94 = 04 — 05, and the additional resonant mode 62 = —Q)^i'p2 — — ^203). 
Employing canonical quantization of this circuit, we find the effective Hamiltonian 

H = Ho+ ^ njjjja]aj + 9jN{a] + hj), (4) 

i=l,2 j=l,2 

describing the inductively shunted junction, Hq, coupled to two resonant modes j = 1, 2 with 
coupling strengths gi and g2, respectively. 

Theory fits to experimental data. Design and fabrication of the fluxonium system only allow 
for imprecise estimates of the system parameters. Thus, the comparison between experimental 
data and theoretical prediction requires the fitting of theory curves to determine the system 
parameters with more accuracy. The parameters at our disposal are: Ej, Eq, and El (for both 
the inductively shunted junction and the extended fluxonium models). In addition, the extended 
fluxonium model takes the resonant mode frequencies and coupling strengths ciJi 2 and (71 2 as 
input. Fits are obtained by extracting the center frequencies from the experimental data and 
employing a least-squares fit algorithm. 

Fit to inductively shunted junction model. A simultaneous fit to the full flux-dependence of 
the 0-1 and 0-2 transitions around the zero-flux point fully determines the fluxonium parame- 
ters Ec, Ej, and E^ (see Table[T]for the obtained parameter values). A comparison between the 
resulting theory prediction of higher transitions can then be used as a consistency check. While 
the agreement for the 0-1 and 0-2 transitions is good, we find systematic deviations for higher 
levels. The reason for these deviations lies in the effect of the additional resonance on the 0-2 
transition: the additional resonance leads to significant frequency shifts of the 0-2 transitions. 
Ignoring this effect leads to a systematic error in the estimation of the fluxonium parameters. 

Fit to extended fiuxonium model. For best agreement, both resonator and additional resonant 
mode are taken into account. Using the full experimental data we obtain a fit for the extended 
fluxonium model, which shows excellent agreement with the data. The resulting parameter 
values are given in Table [T] 
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inductively shunted extended fluxonium model 
junction model 



Eclh 2.39 2.47 

Ej/h 8.93 8.97 

El/H 0.52 0.52 

a;i/27r n.a. 8.18 

a;2/27r n.a. 10.78 

c/i/27r n.a. 0.135 

c/2/27r 0.324 



Table 1: Fluxonium system parameters obtained from least-squares fits to the inductively 
shunted junction and the extended fluxonium model. All values are given in GHz. The coupling 
constants are expressed in terms of the coupling strength for the fluxonium 0-1 transition. 
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